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In this paper we investigate the role played by dark energy in the form of Generalized 
cosmic Chaplygin gas in an accelerating universe described by FRW cosmology. We have 
tried to describe the model from the theoretical point of view of a field, by introducing a scalar 
field <fi and a self interacting potential V ((/)). The corresponding expressions for the field are 
obtained for the given model. Statefinder parameters have been used to characterize the 
dark energy model. Plots have been generated for characterizing different phases of universe 
diagrammatically and a comparative study is performed with the Modified Chaplygin gas 
model. As an outcome of the study, Generalized cosmic Chaplygin gas is identified as a 
much less constrained form of dark energy. 
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I. INTRODUCTION 

At the turn of the last century a new feather was added to Edwin Hubble's discovery of the 
expanding universe, when it was confirmed from the observations of type la supernovae and CMBR 
data that the expansion is having an accelerating nature associated with it [2] - It marked the 
beginning of a new era in cosmology, as it gave rise to the new concept of the mysterious dark 
energy(DE) [3]. DE, as the name suggests is an invisible form of energy, possessing a large negative 
pressure component, with the help of which it is able to bring about the cosmic acceleration. Hence 
this type of matter violates the strong energy condition, i.e., p + 3p < 0. 

Now, soon after this attribution, the entire cosmological society began to investigate for a 
suitable model that can effectively play the role of DE. The simplest form of DE was found to be 
a tiny cosmological constant, A, which gave rise to the simplest model of the accelerating universe, 
the ACDM model. DE associated with a scalar field 4] is called quintessence, and provides an 
effective model of the accelerating universe. As time rolled on various Chaplygin gas models came 
into existence. Extensive research saw Chaplygin gas (CG) |5[ pave way to Generalized Chaplygin 
gas (GCG) jfj and Modified Chaplygin gas (MCG) Q, Q|. The MCG equation of state is given by 

P = A P -^ (1) 

where p and p are respectively the pressure and energy density and < a < 1 , A and B are positive 
constants. 

In 2003, P. F. Gonzalez-Diaz [9] introduced the Generalized cosmic Chaplygin gas (GCCG) 
model. The striking factor of the model, which made it unique, was that it can be made stable 
and free from unphysical behaviours even when the vacuum fluid enters the phantom era. In the 
previous studies related to DE corresponding to phantom era it was seen that Big- Rip was the final 
destiny, since the time derivative of scale-factor blows to infinity in finite time. For the first time 
P. F. Gonzalez-Diaz introducing the GCCG model showed that Big Rip singularity can easily be 
avoided. Hence in such models there is no requirement for evaporation of black hole to zero mass. 

The equation of state of Generalized cosmic Chaplygin gas is, 



P = ~P 



C+{p 1+Q -C} i (2) 



where C = — 1, with A being a constant that can take on both positive and negative values, 
and — L < co < 0, L being a positive definite constant, which can take on values larger than unity. 

As numerous DE models began appearing in the scene, it was of utmost necessity to devise a 
method that would both qualitatively and quantitatively discriminate between various DE models. 



The statefinder parameters [101 ] were introduced in this context by Sahni et al in 2003. These 
parameters are able to discriminate between different DE in a model independent manner. The 
diagnostic pair is constructed by the scale factor a(t) and its derivatives as follows: 

where a is the scale factor of the universe, H is the Hubble parameter, a dot (.) denotes differen- 
tiation with respect to the cosmic time t and q is the deceleration parameter given by, 

In fact for ACDM model this pair take the values r = 1, s = 0. For any DE model the trajectory in 
the r — s plane is obtained and the distance of this trajectory from the point (1, 0) is noted. This 
produces the required discrimination. 

In section 2, all discussions are valid in general for open, closed and flat model of universe, i.e., 
for k = ±1, 0. But in section 3, only the simple case of spatially flat universe (k = 0) is considered. 



Ill ], that the flat model of universe follows 



In fact it is confirmed from various CMB experiments 
spontaneously from the early inflationary phase, thus justifying our consideration. 

In this work we have tried to describe the GCCG model from the theoretical point of view 
of a field, by introducing a scalar field eft and a self interacting potential V((ft) with the effective 
Lagrangian 

L = \ft- V(<ft) (5) 



In 



it is shown that the flat Friedmann model with Chaplygin gas can be described equivalently 



in terms of a homogeneous minimally coupled scalar field (ft, by fitting the FRW equations for 



Chaplygin gas into Barrow's scheme 



121 ] . The expressions for homogeneous scalar field ( ft(t) and 



potential V((ft) describing the Chaplygin cosmology was obtained by Kamenshchik et al 
In the work was repeated considering MCG as the DE model. In this work we have obtained the 
corresponding expressions for GCCG model as a generalization of the previous works on Chaplygin 
gas cosmology. 

The paper is organized as follows: In section 2 the expressions for <ft(t) and V{<ft) are obtained in 
terms of the scale factor a(t) and the corresponding plots are generated to show how V{(ft) varies as 
the scale factor varies. In section 3, the evolution of the universe is studied in the {r, s} plane for 
the entire physically realistic history of the universe and compared to that of the other Chaplygin 
gas models, especially MCG. In section 4, a complete graphical analysis is performed. In section 5 
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a comparative study between MCG and GCCG is performed and finally the paper ends with some 
concluding remarks in section 6. 

II. GENERALIZED COSMIC CHAPLYGIN GAS IN FRW COSMOLOGY 

The metric of a homogeneous and isotropic universe in the FRW model is 



ds 2 = dt 2 -a 2 {t) 



dr 2 



1 — kr 2 



+ r 2 (dO 2 + sin 2 



where a(t) is the scale factor and k(= 0, ±1) is the curvature scalar. 
The Einstein field equations are 

» 2 k 1 



a 2 a 2 3 



P 



and 



a 1 , 

- = --(p + S P ) 

a o 



(6) 



(7) 



(8) 



where p and p are the energy density and the isotropic pressure, respectively. Here we choose 
8vrG = c = 1. 

The energy conservation equation is 



p + 3-(p + p) = 



(9) 



Using equation ([2]) in equation Q, we have the solution of p as 



P 



C+{\ + 



B 



3(1+q)(1+<j) 



a 



i 



1 

1+a 



(10) 



where B is an arbitrary integration constant. 

Using eqns. (j71 [HI and [TUj) we get the following differential equation for k = and a = 1. 



a(t) 



o(t) 
6 



a(t) 
a(t) 



9 



o(0 
a(t) 



C+ 9 



a(t) 
a(t) 



C 



(11) 



This being a highly non-linear differential equation, it is very difficult to find an explicit expres- 
sion of the scale factor by the known mathematical methods. 
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Now for small values of scale factor a(t), we have, 



B (l+w)(l+a) 

< 12 > 

which is a very large value and corresponds to the universe dominated by an equation of state, 
p = —\, where Z = a (1 + uj) + oj is a constant. 
For large values of scale factor a(t), we have, 

+ (13) 

For this value of p, the value of p is found to be, 

p ~ - (C + l) 1 ^ - _p (14) 



which corresponds to an otherwise empty universe with a cosmological constant (C + 1) 1 + a . 

We know that an accelerating universe should have a negative deceleration parameter, q. So it is 
clear from the expression of deceleration parameter (q = — -^jp:) that, a > 0. Some straightforward 
algebra using equations (|8l and [TOl) . gives, 



i 



X [l-2CX~—) < 3 (15) 



where X = 1 H — 3(1+1 f )(1+Q) • Here the transition from decelerating to the accelerating universe 



occurs at a = ^ — ^ — J . Considering the sub- leading terms in eqn. (llOh we obtain for large values 
of scale factor, the following expressions for energy density and pressure: 



p * + y c) (i6) 

C~ (1 + a) 

and 



-3- ( y - C) 

l+a -| i- 1 

C 1 ^ (1 + a) 



(17) 



where Y = C + {I + 



r} 



a 3(l + ^)(l + a) 

Now we consider this energy density and pressure corresponding to a scalar field (p having a 
self-interacting potential V(4>)- the Lagrangian of the scalar field is given by, 



H = ~ v(<t>) 



(18) 



6 



The expressions for energy density p§ and pressure p^ for the scalar field are as follows, 



B 



3(14«)(l+w) 



i 

l + LJ 



1 

1 + Q 



and 



(19) 



1 ■ 



p<p = -J A - v{4>) 



C+\l + 



B 



3(l+a)(l+o;) 



C+ 1 + 



B 



3(H-a)(l+w) 



1+u 



= -Y-ttz [c + (Y -cy 
For a flat universe, i.e., k = 0, we have, 



(20) 



B 



3(l+a)(l+w) 



l 

1+w 



1- <C+il+ 



B 



j3(l+a)(l+w) 



1+w 



and 
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3(l+o)(l+w) 













i 




f . 





y{c + (y-c)^} 



(21) 



v(<t>) = - 



3(l+a)(l+a;) 



l -| 

1+w 



1 

1+a 



a 3(l+a)(l+w) 



C+U + 



5 



; 3(l+a)(l+aj) 







i 






y l+a 


n 


i. 


2 



1 + 



y{c + (y-cr w } 



(22) 



From the equation (|2ip . we have the following relation between scale factor a(t) and scalar field 



(-1 
C + (1 + Bav~ 3 ( 1 + w )( 1 + a ))~ 
1 — 
C+ (1 + B au -3(l+u ) )(l+a))i+ 



(23) 




where C is the integration constant. 
Now 2V(cf>) = p-p = p- C+(pl+ °" C) 



Case I: 

When p — > oo, i.e., a — > 0, — > oo 
Case II: 

Whena^oo, V{4>) ->• (C + l) 1 ^ 

For example, 

When C = 1 and a = 1, V(0) -»• \/2 as a — > oo 
When C = -1 V(0) ^ as a -»• oo 

When C = V(#) -»• 1 as a ^ oo 



Description of figures 

Figs.l : The potential is plotted against time for different values of a. The values of a used 
are 0.1 and 0.8. The value of uj is taken as —0.5. 



Fig. 5 : The potential is plotted against time for different values of a. The values of a used are 
0.1 and 0.8. The value of cv is taken as —1.5. 
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III. GEOMETRICAL DIAGNOSTIC USING STATEFINDER PARAMETERS 



The statefinder parameters are defined as follows, 



r "dp S 55 3 J- 1/2) (24) 



where H is the Hubble parameter and q = — is the deceleration parameter. Trajectories in the 
r — s plane characterize different cosmological models. The ACDM model corresponds to the fixed 
point {1, 0} in the r — s plane. The deviations of the trajectories of different cosmological models 
from this point in the r — s plane gives an effective comparison of the model from the standard 
model. So we proceed towards the statefinder analysis of the present model, GCCG. 

For Friedmann model with a flat universe we have 



and 



aH 2 2 2p 

From the equation (|24p . we have, 



l + l 1 (26) 



-= 1+ K 1+ ?)0 • -=( 1+ i)0 (27) 



We get the ratio between p and p as, 

p 2(r-l) 



p 9s 

For GCCG the velocity of sound can be given as, 



(28) 



VS ~ dp p + (pl+a _ C ) w+1 (29) 

From equations (|27]h (|28p and (|29p . we get the following relation between r and s 

18s 2 (9a) « (r - l)+(9s) 1+ " 2a (r - 1)+4q (9s) « (r - l) 2 = w (1 + a) (2 - 2r - 9C,s) 1+ - (9s + 2r - 2) 

(30) 

Now we consider three different cases: 
Case 1: When uj = — | 



For this case equation (j30]l becomes 

18s 2 (r - l) + 18a (r - l) s + 4a (r - if (2 - 2r - 90s) + 81s^ (1 + a) (9s + 2r - 2) = (31) 
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From the above equation it is quite clear that there is only one asymptote parallel to the s-axis, 



A x = 1 - 18C + 2a - l8Ca + 4C 2 a + a 2 

B x = (-17C 2 + C 2 a) 

Cx = a + C 2 a + a 2 + C 2 a 2 

Dx = -2 + 108C - 1944C 2 + 11664C 3 - 12a + 540Ca - 8412C 2 a + 46872C 3 a - 39600C 4 a + 
11016C 5 a - 30a 2 + 1080Ca 2 - 13560C 2 a 2 + 58968C 3 a 2 - 78648C 4 a 2 + 23112C 5 a 2 - 2448C 6 a 2 - 
40a 3 + 1080Ca 3 - 9648C 2 a 3 + 23976C 3 a 3 - 35976C* 4 a 3 + 11448C 5 a 3 - 2432C 6 a 3 - 30a 4 + 540Ca 4 - 
2544C 2 a 4 + 216C 3 a 4 + 3000C 4 a 4 - 648C 5 a 4 + 144C 6 a 4 - 12a 5 + 108Ca 5 + 12C 2 a 5 - 72C 4 a 5 - 2a 6 

Case 2: When u = — 1 

For this case equation (i30"j) reads 





where 



2(r-l) s + a + 



2a (r - 1) 
9s 



+ (1 + a) (9s + 2r - 2) = 



(32) 



Here the only asymptote parallel to the s-axis is r = 1 



9(l+a) 



. Here the asymptote intersects the 




Case 3: When to 



= -2 



Here the equation (]30p transform into 



9(1 + a) 2 s (2 - 2r - 9Cs) (9s + 2r - 2) 2 + 9s 2 (r - 1) + 9a (r - 1) s + 2a (r - 1) 



(33) 
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where, 

A 2 = 8VC + 4C 3 / 2 + 16y/Ca - 4Ca + 8C 3 / 2 a + 7VCa 2 - ACa 2 + 4C7 3 / 2 a 2 
B 2 = 1 + 2C + a - VCa + 2Ca 
C 2 = 2C + 6Ca + 5Ca 2 + Co? 

D 2 = 128C 3 / 2 + 1344C 5 / 2 + 384C 7 / 2 - 128C 9 / 2 + 768C 3 / 2 a - l92C 2 a + 8064C 5 / 2 a - 192C 3 a + 
2304C 7 / 2 a + 384C 4 a-768C 9 / 2 a+2016C 3 / 2 a 2 -960C 2 a 2 + 20352C 5 / 2 a 2 -960C 3 a 2 +6624C' 7 / 2 a 2 + 
1920C 4 a 2 - 1920C 9 / 2 a 2 + 2944C 3 / 2 a 3 - 2160C 2 a 3 + 27648C 5 / 2 a 3 - 2848C 3 a 3 + 11136C 7 / 2 a 3 + 
3840C 4 a 3 - 2560C 9 / 2 a 3 + 2520C 3 / 2 a 4 - 2640C 2 a 4 + 21672C 5 / 2 a 4 - 4704C 3 a 4 + 10944C 7 / 2 a 4 + 
3840C 4 a 4 - 1920C 9 / 2 a 4 + 1200C 3 / 2 a 5 - 1728(7 2 a 5 + 9552C 5 / 2 a 5 - 3744C 3 a 5 + 5760C 7 / 2 a 5 + 
1920C 4 a 5 - 768C 9 / 2 a 5 + 250C 3 / 2 a 6 -480C 2 a 6 + 1896C 5 / 2 a 6 - 1120C 3 a 6 + 1248C 7 / 2 a 6 +384C 4 a 6 - 
128C 9 / 2 a 6 
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Description of figures 

Fig. 3 : The statefinder parameters are plotted against each other for different values of a. The 
values of a used are 0.5, 1 and 5. The value of C is taken as 2. 

Fig.4 : The statefinder parameters are plotted against each other for different values of C. The 
values of C used are 6 and 20. The value of a is taken as 0.5. 

Fig. 5 : The statefinder parameters are plotted against each other for different values of C and 
a. The values are furnished in the figure. 

IV. GRAPHICAL ANALYSIS 

In the figures 1 and 2, the scalar field potential V{4>) is plotted against time for different values 
of the parameters involved. Two different plots are generated for different values of to. It is evident 
that the scalar field potential decreases as the scale factor increases. In an accelerating universe the 
scale factor increases gradually with time. So it is clear that the potential decreases with 
time. Plots with different values of a shows that the potential assumes larger values for smaller 
a. On the contrary a comparison between figs. 1 and 2 indicate that with a decrease in the 
value of oj, there is corresponding decrease in the value of potential. 

In figure 3, the statefinder parameters r and s are plotted against each other for three different 
values of a. It is seen that for r = 1, all the three curves tend towards s = 0, thus producing a 
correspondence with the ACDM model. With an increase in the value of a, there is a gradual shift 
of the trajectories towards the right hand side of the plot, thus showing an increasing r tendency. 
In this plot C = 2 is kept constant throughout. 

In Figure 4, the r-s trajectories are obtained for two different values of C, keeping a fixed. It 
is seen that for C = 6 (red curve), s = more or less corresponds with r = 1, whereas for C = 20 
(blue curve), s = does not exactly correspond to r = 1, thus showing a greater deviation. Thus 
it can be concluded that as the value of C increases, there is greater deviation of the 
GCCG model from the standard ACDM model. Since C can take both positive and negative 
values, r-s plots for negative values of C should be obtained and studied. This is accomplished in 
fig. 5. Three different values of a are considered and the plots are generated. 
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V. A COMPARATIVE STUDY BETWEEN THE ROLE OF GCCG AND MCG IN 

ACCELERATING UNIVERSE 

In 8J] the authors studied the role of MCG in the accelerating universe. In this section we 
undergo a comparative study between this study and the study in |8||. We revisit the plots for the 
scalar field potential and statefinder parameters for MCG as witnessed in [8J] in the figs. 6, 7, 8, 9 
and 10. These plots when compared to the corresponding plots of GCCG exhibits a few similarities 
and dissimilarities. We study them separately below. 

A. Similarities 

We see that the plots 1 and 2 for GCCG almost show similar characteristics as the plots 6 and 
7 for MCG, taking A = 1/3. In these figures the potential decrease with increase in scale factor. 

B. Differences 

1) In the plots 8 and 9 for MCG, taking A = 1, the curve behave quite differently from the plots 
1 and 2 for GCCG. In case of MCG the potential increases steeply for a while before becoming 
stagnant at a particular value of V. But in case of GCCG there is a constant decrease in potential 
with increase in a, as seen before. 

2) Fig 10 shows the r — s plot for MCG. When compared with the r — s plots for GCCG (figs. 
3, 4 and 5), clear differences in trajectories are observed. This helps us to discriminate between 
the two DE models. 

3) In fig. 10, the portion of the curve on the positive side of s, which is physically admissible is 
only the values of r greater than {1 + |^4 (1+^4)} (for A = 1/3, it becomes equal to 3). Incidentally 
the portion of the curve between r = 1 and r = {1 + §4_(1 + ^0} on the positive side of s is not 
admissible and that portion of the curve is not shown in 
restricted region constant B becomes negative, which is inadmissible. Therefore the curve on the 
positive side of s starts from the radiation era and goes asymptotically to the dust model, but the 
portion on the negative side of s shows the evolution right from the dust state (s = — oo) to the 
ACDM model (s = 0). Hence the negative and the positive region complement each other, and the 
total curve (both positive s and negative s portions) represents the evolution of the universe from 
the radiation era to the ACDM model. 



8]. The reason behind this is that in this 
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Variation of V against a 



a=l 
A=l/3 



MCG case 



Variation of V against a 



a=0.6 
A=l/3 



MCG case 



Fig. 6 

Variation of V against a 



a=l 
A=l 



MCG case 



Fig. 7 

Variation of V against a 



a=0.6 
A=l 



MCG case 



Fig. 8 



Variation of s against r 



A=l/3 / 
a=l 1 




MCG case 





Fig. 9 



Fig. 10 



But in case of GCCG we see that there is no such restriction on any of the parameters C, u 
and a, since they do not have any restriction on either magnitude or sign. So the entire portion of 
the curves in figs. 3, 4 and 5 is physically admissible and there is no restricted portion, as it is in 
case of MCG. This is a striking difference! In fact it makes GCCG a much less constrained 
form of dark energy. 
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VI. CONCLUSION 

In this work we have considered Generalized cosmic Chaplygin gas and tried to find the role 
played by it in an accelerating universe described by FRW cosmology. To accomplish this we have 
compared the GCCG model with a scalar field, and determined its potential. Plots were generated 
for the determined potential. Statefinder parameters were studied in detail in order to discriminate 
the DE model from other DE models. We have undergone a comparative study between the GCCG 
and MCG models as far as their role in the accelerating universe is concerned. The trajectories 
of the r — s plots helps us to discriminate between the two models. For A = 1, it was found that 
MCG shows an increase in the scalar potential in a certain region, as a increases. But in case of 
GCCG, irrespective of whatever value is assigned to the parameters, the potential always show a 
downhill progress. It was witnessed that unlike MCG there is no restricted region in the r — s 
plot for GCCG, thus making every portion of the curve physically admissible. Finally it should be 
stated that GCCG describes the universe from radiation era to ACDM era without constraining 
its parameters considerably. This makes GCCG a much less constrained form of Dark 
energy. 
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